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The behavior of a family of mean-field glass models is reviewed. The models are analyzed 
by means of a Langevin-based approach to the dynamics and a Replica theory computation 
of the thermodynamics. We focus on the phase diagram of a particular model case, where 
glass-to-glass transitions occur between phases with a different number of characteristic time- 
scales for the relaxation processes. The appearence of Johari-Goldstein processes as collective 
reorganizations of sets of fast processes is discussed. 



Theoretical modeling of glassy systems is a widespread topic. Different important 
theories have been introduced along the years like, e.g., the "free volume" [1], 
the "entropic" [2] and the "random first order" [3] theories, to mention a few. 
Nevertheless, a comprehesive theory both being analytically treatable and yielding 
reliable quantitative predictions above and below the glass transition has yet to be 
devised and appears to be a very complicated and challenging aim. In front of such 
limitations a theoretical approach based on the mean-field approximation, where 
statistical fluctuations of microscopic observables are neglected, helps pointing out 
a way to enforce realistic approaches and indentifying physically relevant concepts. 

We present the study of a family of mean-field models with time- independent, 
i.e., quenched, disorder. The typical feature of amorphous systems is the impossi- 
bility to reach states at the lowest feasible energy, thus preventing crystallization. 
The system, undergoes some kind of "frustration" : because of dynamic arrest the 
global set of the energetic contributions due to the interactions among the glass 
former constituents cannot be simultaneously minimized. In the present case the 
frustration is a direct consequence of the quenched disorder. More generally it is 
self-generated by the inner geometry of the material and/or by the complicated ex- 
change of interactions. The quenched disorder is not a necessary ingredient Q then, 
even though it makes the problem more easily tractable. The motivations for this 
study are manifolds. First, applying the replica method and using the concept of 
Replica Symmetry Breaking (RSB), the analysis of thermodynamic and dynamic 
properties can be carried out analytically. Then, it is possible to develop and check a 
multi-timescales equilibrium dynamics consistent at all times, including the asymp- 
totic limit. Further, the model displays a very rich phase diagram with different 
glass and spin-glass phases, allowing for a theoretical analysis of the phenomenon 
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Tabic 1. Phases and order parameters for models with quenched disorder. 



PHASE 


OVERLAP DISTRIBUTION 


Paramagnet /Fluid 


P{q) = 


%) 






Glass 


P(q) = mS{q - q ) H 


-(1- 


m)S(q — 


91 ) 


Spin-Glass* 


P(q) = w 5(q - qo) + 


Hl) 


+ w\5(q 


-91) 



* The function P is continuous on the support ]qo : qi[, wo and w± are weights of the <5's (wq + wi < 1). 



of polyamorphismj^] Eventually, the identification of RSB's in the thermodynamics 
with time-scale separations in the dynamics [TT] provides an useful theoretical tool 
to study the interrelation between primary (a), secondary (or Johari-Goldstein, 
/3jg) and tertiary processes (/3f ast ,7). In this paper we will deepen the last two 
aspects and their possible implications for real structural glasses. 

We first briefly introduce the model and sketch the computation of its thermody- 
namics within the Replica theory, emphasizing the nature of the order parameter 
and its change in behavior accross qualitatively different amorphous phases. The 
model Hamiltonian is: 

ii<...<« s ii<...<i p 

where J® it (t = s,p) are uncorrelated, zero mean, Gaussian variables of variance 
J t 2 t!/(2A^* _1 ) and crj are N "spherical spins" obeying the constraint J^of = N. 

In a complex Free Energy Landscape (FEL), such as the one representing an 
amorphous system, the numerous valleys, i.e., the "glass states", can be more 
or less correlated among them and a hierarchy can be established based of their 
relative correlationj^Denoting by (. . .) a the thermal average over the configurations 
belonging to state "a", the following overlap order parameter is defined as the 
correlation between two states (a and b): 



Qab= j^J2{vt)a(cri) b (2) 
i=i 

To be precise, the complete order parameter is the probability distribution P(q) of 
the values of q p3] . Depending on the shape of P(q) one can identify a specific phase 
of the amorphous system. In table [T] we summarize the most common behaviors 
known in literature. 

The replica theory for mean-field disordered systems is applied to compute the 
free energy functional |12| [13] : 

1 1 1 1,n 

- p<S> = -(1 + In 2tt) + - lim - V g{q ab ) + In det q (3) 

I l n^O n * — ' 

ab 

where q = {q a b} is the Parisi overlap matrix. The model is specified by the function 
9{q) = q s l 1 s/s + q p Vp/p, with fj, s = sf5 2 J^/2. For a generic RSB Ansatz with R 
breakings the elements of the Parisi matrix take values = go < Ql < • • • < 1R < 



2 Many examples of polyamorphism are available in nature (and in literature). For example, the change in 
the kinetics of the coordination between molecules, occurring in vitreous Germania and Silica [5] or the 
sharp density change taking place in porous silicon [B], as well as in undercooled water 0. Very recently 
polyamorphism in Ethanol 0, Laponite [9] and star polymer mixtures 1101 has been observed. 
J A well known example is the symbolic dynamics through the Potential Energy Landascape, where intra- 
basin processes have a high correlation and inter-basin processes have a low correlation |16l 1171 . 
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qn+x = 1 with relative multiplicities n = too > toi > . . . > ran > mR+l = 1- As, in 
the Replica computation, n — > 0, the parameters m r acquire real values (€ [0, 1]) 
[T3] and one can express the set of q and m values as a (step) function q(x). 
Here we are interested in structural glass. We will, thus, take into account model 
cases displaying phases with one and two step RSB, whose overlap functions are 
schematically represented on the left hand side of Fig. [T] These can be qualitatively 
connected with real glass formers in which only primary (R = 1) or also secondary 
(R = 2) processes are present. [^]Such glass models are realized taking s > 2 and 
large p — s [201 EI] • Eq. (|3j for a R = 2 RSB phase can be written as 

- 2/3$ = 1 + In 2vr + g(l) + m 2 [g(q 2 ) - g{q 1 )} + mi[#(gi) - g(q )} (4) 

+ In X2 In In 1 (5) 

m 2 Xi rni xo Xo 

with X2 = 1 - 92, Xi = X2 + m 2 {q2 - qi) and X o = Xi + ™i(<7i ~ qo)- 
In the s + p models there are also different solutions (depending on the values of 

s and p and of T, J s and J p ), displaying both continuous and discontinuous (and 

mixed) overlap functions. |1 2|, I13j ^] 

The order parameter function (i.e., the set of values of m's and g's) is obtained 

by solving the following set of self-consistency equations: 

5(92) - g{qi) = (92 - 9i) U(9i) — j - ^2 ln — (6) 

L m 2 xii 171$ xi 

g{qi) - 5(90) = (91 - 90) |a(<? ) — j - \ In — (7) 

L miXoJ rnf xo 

A(qo) = % A( qi )-A(q ) = ^^; A(q 2 ) - A( Q1 ) = ^± (8) 
Xo X0X1 X1X2 

with A(q) = dg{q)/dq. The thermodynamics of the 1RSB solution is obtained from 
the above equations setting q 2 = q%. 

The value of the overlap corresponds to a given correlation among states, cf. Eq. 
([2]). The three levels function displayed in the 2RSB solution corresponds to a pre- 
cise hierarchy in the organization of the thermodynamically relevant glassy states, 
consisting in groups of states (clusters) and groups of gropus of states ("meta"- 
clusters). Two states whose overlap is q 2 belong to the same cluster. Two states 
whose overlap is q± do not belong to the same cluster but to the same meta-cluster. 
Eventually, two states whose overlap is qo (usually equal to zero in absence of ex- 
ternal forces of fields) belong to different meat-clusters. 

The dynamics of the model is Langevin. Using a Martin-Siggia-Rose path-integral 
formalism one can reduce the equations of motion to a single variable formulation. 
Details can be found in Ref. [21]. The most important two-time observables are 
the correlation and response function 



C(t,t w ) = (a(t)a(t')h G{t ^ tw)=8 ^) f - tw (9) 

Let us take an amorphous phase with a generic number R (= 1, 2, . . .) of time-scale 



2 From the point of view of Replica calculation we stress that a thermodynamically consistent example of 
a 2RSB phase has not been realized in models other than the s + p spherical models 1131 . 
1 Usually, a glass phase is associated with discontinuous steps in the overlap, corresponding to a sharp 
separation of time-scales. A spin-glass phase is, instead, characterized by a fully continuous q(x). 
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Figure 1. L.h.s.: order parameter step function q(x) for the 1 and 2 RSB thermodynamic glassy phases. 
R.h.s.: Correlator vs. time for systems with processes relaxing on one and two well separated time-scales. 
A correspondence between the properties of the overlap function in the static solution and the behavior of 
the correlation function in the dynamics is shown (see also Eq. | |12| l). 

bifurcations and assume that equilibrium is obtained in each (completely) disjoint 
time-sector, i.e. t w — > — oo. Using, e.g., a multiple scale analysis one assumes that 
the correlation function C(i), as well as the response G(t), can be represented as 
the sum of R + 1 distinct terms each depending on a time variable, to <C ■ ■ • T r <C 
... <C tr, describing the motion in a given time-sector: 



R 



C(t) = Y, C r (l 



(10) 



Considering a time-sector a means to probe the dynamics on times t ~ r a . We can 
split off the a-sector function C a taking the ordered limit 



lim^oo = lim . . . lim , 



T ^OO 



(11) 



with the prescription r a jt = O(l), r r<a /t — > and T r>a /t — > oo. In practice, all 
contributions C r to the correlation function with r < a correspond to processes 
already thermalized at the observation time t, whereas all contributions whose 
index is larger than a represent processes that are frozen at time t. The interesting 
processes under probe are those relaxing on characteristic times r a ~ t. In the 
above formulation, the asymptotic value of the correlator is, then, C{t) — > q r , and 
we have the condition: 



lim t . 



r-l 



^C s (r s ) + ^C s (r s ) 



Vr = 0, 



,R 



(12) 



A schematic behavior of C{t) for the cases of our interest, R = 1,2, is plotted 
on the r.h.s. of Fig. [T] next to their overlap counterparts. In the top part we 
have the thermodynamic order parameter (step) function q{x) displaying a single 
discontinuity at x = mi. The two segments of the step function (qo,qi) can be 
linked, in the dynamics, to the two plateaus of the relaxation function in ordinary 
glass formers, in cases where secondary processes play no role. In the bottom part 
the 2RSB case is sketched, i.e., two discontinuities in q{x) (thermodynamics), or 
two time-scale separations in C (dynamics). This is likely to be the mean-field 
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reduction of a glass with secondary processes. 
The response function on multiple separated time-scales reads 

R 
r=0 

where each function G r varies only in the corresponding sector r, r r ~ 0(t) and 
vanishes in all sectors with s < r. The function G r represents the response of the 
system to a perturbation in the time sector labeled by r, i.e., the response due to 
all degrees of freedom which have not equilibrated in previous sectors. 

Working with the Fourier transforms of the correlation and response functions 
and defining the kinetic coefficient r~ 1 (w) = idujG^ 1 (lo) , the dynamical stability 
is guaranteed by the requirements T(u r ) = 0, \/r = 1, . . . , R and T(luq) > 0, as the 
ordered limit lim^^o • • • hm^ ^o is performed. In our R = 2 case the conditions 
can be written as 

A'fe) = 1/x!, A'(9i) = l/xl A! (go) > 1/Xo ( 14 ) 

and provide the equations for the asymptotic dynamic solution. We stress that the 
solution to Eq. (14) do not coincides with the static solution, Eqs. ([6])-(|8]). This 



is typical of systems undergoing a dynamic arrest before reaching a temperature 
where they can undergo a thermodynamic phase transition. The thermodynamic 
transition in these spin-glass inspired mean-field models for the glass is, instead, 
the Kauzmann transition (at Tr), whereas the dynamic transition (at Tj) is equiv- 
alent to the dynamic arrest transition predicted, e.g., in Mode Coupling Theory 
(MCT). In real experiments it corresponds to the crossover temperature at which 
the separation of time-scales of slow and fast processes accelerates. Q 

A straightforward link with schematic models in MCT [18] can be drawn, start- 
ing from the observation that the dynamic equations in random spherical models 
are equivalent to the MCT equations [19] at high temperature, where time trans- 
lational invariance (TTI) holds and G and C are connected by the fluctuation- 
dissipation theorem (FDT), G{t — t') = —f39{t — t')dtC(t — t'). In our model case, 
thus, if we take a memory kernel depending on the correlator <fi (in MCT notation) 
as m((f)) = /u s s_1 + n P (f) p ^ 1 = A(0), the mode coupling equations describe the 
Langevin dynamics of the s + p model and the overlap is identified with the non- 
ergodicity parameter: q = lim^oo </>(*)• The two dynamics differ, instead, below 
Td, since the global TTI breaking is implicit in the random model dynamics, and 
FDT does not apply anymore in the above form abovej^] 

In mean- field models, unlike real glasses, the configurational entropy S c = log J\fj 
is a true state function. It can be formally computed as the Legendre transform of 
the total free energy <£, with / and (3m as conjugated variables: 

SJf; T) IN = min [-/3m$(m; T) - (3m f] (15) 
m L 

where m = trr- The configurational entropy of the thermodynamic solution is 
subsextensive (S c (f eq )/N — > 0, ./V — > oo, Kauzmann point). By maximizing S c vs. 



J We notice that the experimental, calorimetric, glass temperature T g is not defined in mean-field systems. 
Indeed, this is a property connected with the falling out of equilibrium of activated processes (hopping 
among valleys), whereas in mean-field metastable states are surrounded by infinite barriers (as N — » oo). 
T g lies, undetermined, between Tx and T^. 

2 For details on the generalization of equilibrium dynamics in the solid amorphous phase see Ref. |21|V 
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Figure 2. Phase diagram of the 3 + 16 model with both 
dymanic (dotted) and thermodynamic (full) transition 
lines. Four phases are present: paramagnetic (PM) at 
high T and three glass phases, here termed 1RSB/, 
lRSB /7 and 2RSB (see text). 



Figure 3. Interpretation of the GGT's in 
terms of metastable states organization. In one 
case a GGT occurs when fast processes slow 
down (and become secondary) and even faster 
processes appear (left, "states fragmentation"); 
in the other case secondary processes show up 
as new, intermediate, processes, between /3f as t 
and a processes. 



/, instead, we find the same values of q r and m r that solve the dynamic equations. 
As a consequence, the dynamic arrest temperature can also be identified by looking 
at the temperature at which an extensive configurational entropy arise. 

In Fig. [2] we show a detail of the (T/ J s , J p / J s ) phase diagram of the (s,p) = 
(3, 16) model around the tricritical point. Both dynamic and thermodynamic (i.e., 
Kauzmann) transition lines are plottedj^] One can observe that descreasing T the 
dynamic transition between glass phases of different nature always precedes the 
thermodynamic one. The dynamic transition temperature Td{J p ) is the highest T 
at which the lifetime of high-lying local states becomes infinite and their number 
grows like exp(S' c (T)) with the size N. The Kauzmann temperature Tk(J p ) is, 
instead, the highest T at which S c (f eq )/N of the global gZass^ minima goes to 
zero. The overlap order parameter jumps from to qea while the free energy $ is 
continuous: ®n q (T K ) = $ g i ass (T^). 

At low T, two glass to glass transitions (GGT) occur, at the lower and higher 
J p . We notice thay they are not exactly of the same kind, in terms of local states 
hierarchy change. We schematically report in Fig. [3] how the states in the 1RSB 
phase reorganize as the system transforms into a 2RSB glass in the two cases. In 
the transition at low J p a local state fragments into a cluster of new local states, 
whereas accross the transition at high J p subsets of uncorrelated local states group 
together in correlated clusters. 

In conclusion, we have examined a mean-field model displaying, in particular, a 
phase whose thermodynamics is described by a 2RSB solution (q(x) = qo,Qi,Q2)- 
Recalling both the thermodynamic and the dynamic properties of this specific 
phase (including the GGT's from other glassy phases) and exploiting the equiva- 
lence between RSB's and time scale bifurcations, we argue that 

(1) changes connecting two local states in the same cluster of states {q = 52) 
are /3f ast (else called 7) processes, 

(2) changes connecting two local states in two different clusters belonging to 
the same cluster of clusters (q = qi) correspond to JG processes, 

(3) changes connecting two uncorrelated states (q = qo « 0) contribute to the 
a relaxation. 



lr The thermodynamic transition, termed "Kauzmann" in the figure, is a so-called "random first order 
transition", with no latent heat but a discontinuous order parameter. This is an example of the mean-field 
scenario behind the mosaic theory [3] 

2 At higher temperature the global minima is a liquid/paramagnetic state. 
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The hierarchical nesting implicit in the present approach hints that fast processes 
have a relevant influence on slow processes, even though taking place on well sepa- 
rated time-scales. This very heuristic observation naturally stimulates a comparison 
with Ngai's Coupling Model (see, e.g., Ref. [23 and references therein). A study 
in this direction is in progress. 
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